Orbital and valley state spectra of a few-electron silicon quantum dot 
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Understanding interactions between orbital and valley quantum states in silicon nanodevices is 
crucial in assessing the prospects of spin-based qubits. We study the energy spectra of a few-electron 
silicon metal-oxide-semiconductor quantum dot using dynamic charge sensing and pulsed-voltage 
spectroscopy. The occupancy of the quantum dot is probed down to the single electron level using 
a nearby single-electron transistor as a charge sensor. The energy of the first orbital excited state is 
found to decrease rapidly as the electron occupancy increases from = 1 to 4. By monitoring the 
sequential spin filling of the dot we extract a valley splitting of ^230 /xeV, irrespective of electron 
number. This indicates that favourable conditions for qubit operation are in place in the few-electron 
regime. 



I. INTRODUCTION 



II. EXPERIMENTAL METHODS 



Silicon-based quantum dots (QD) are recognized as a 
promising system for the implementation of solid state 
quantum computing^"— . While such devices are appeal- 
ing because they can be engineered to show long co- 
herence times^~— , there also exist complications for the 
manipulation of quantum states which arise from the 
multi-valley conduction band present in bulk silicon^-. 
In fact, valley degeneracy is typically lifted in nanode- 
vices due to interfacial stress/strain and electrostatic 
confinement^^^^^. Hence, understanding valley physics 
and its perturbing effects on spin relaxation process has 
become paramount for the practical realization of spin- 
based quantum bits (qubits). In order to assess whether 
a quantum device shows valley splitting compatible with 
use in spin-based quantum computing, orbital and valley 
spectra must be extracted. In particular, it is of inter- 
est to probe the relative magnitude of the valley and 
orbital energy spacing to determine the degree of mixing 
of these states^^ and ultimately identify the most appro- 
priate strategy for qubit operation^^. 

In this work, we use pulsed-voltage spectroscopy to in- 
vestigate the excitation energy spectra of a nearly-closed 
silicon metal-oxide-semiconductor (MOS) QD in the few- 
electron regime {N < 4). In a magnetic field we observe 
Zeeman shifts of the ground state, valley excited state, 
and first orbital excited state that allow us to determine 
the sequential spin filling of the A^-electron states. Inte- 
restingly, the energy of the first orbital excited state is 
seen to decrease rapidly as the dot occupancy increases, 
whereas a valley splitting of ~230 jieV is found not to 
depend on electron number. This has the significance 
of guaranteeing the existence of a well-defined spin- 1/2 
qubit Hilbert space, sufficiently separated from higher- 
energy excitations. 



Fig. [TJa) shows a scanning electron microscope (SEM) 
image of a device identical to the one under study; in 
the experiments presented here, only the (false) coloured 
gates are used. The device works as a single-lead quan- 
tum dot (upper structure) integrated with a single- 
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FIG. 1: (Color online) (a) SEM image of a quantum dot inte- 
grated with a SET and the measurement set-up. Only colored 
gates are used in the experiments. (b)(c) Schematic energy 
diagrams during pulsing of gate P with two different Vp off- 
sets. 
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electron transistor sensor (lower structure). It was fab- 
ricated using multi-layer AI-AI2O3-AI gate stack tech- 
nology and the detailed fabrication process has been de- 
scribed in Refs.^'^^. A quantum dot is formed between 
the left and right barrier gates (LB, RB) and is inde- 
pendently controlled by the plunger gate (P). By apply- 
ing a positive voltage on the lead gate (LD), an electron 
accumulation reservoir is induced at the Si/Si02 inter- 
face. Electrons from the reservoir can then tunnel on 
and off the quantum dot through the left barrier. Next to 
this single-lead quantum dot is a single-electron transis- 
tor (SET) sensor (S) which is used to detect the electron 
tunneling events to/from the quantum dot. The SET top 
gate (ST) forms an electron channel from source to drain 
and the two underlying barrier gates (SLB, SRB) create 
tunnel barriers, forming an electron island in between. 
The measurements were performed in a dilution refrige- 
rator with base temperature of ~50 mK. 

In order to characterize the excitation spectrum^^ of 
the QD, we make use of the technique developed by Elz- 
erman and co-worker^ which combines charge detection 
and gate pulsing. We apply an ac excitation voltage of 
200 /iV at /s=173 Hz to the drain of the SET and mon- 
itor the sensor current Is locked-in to fs at the source 
through a low- noise room temperature preamplifier. Si- 
multaneously, a train of voltage pulses with amplitude 
^puise is applied to gate P, in addition to its dc voltage 
Fp, shifting the energy levels of the dot up and down. 
This pulse train modulates the sensor current via a cross- 
capacitance at a frequency /pulse and the resulting cur- 
rent /pulse at this frequency is measured with a second 
(Pulse) lock-in amplifier, as shown in Fig. [IJa). In or- 
der to maximize charge sensitivity in the detector, we 
employ dynamic compensation, as described in Yang et 
al}^ . We first tune the SET so that the sensor lock-in 
current /s is at the edge of a Coulomb peak, where the 
transconductance dls/dV^T is high (2 nS). We then mon- 
itor the charge state of the quantum dot using the sensor 
signal /s while the SET gate voltage V^t is dynamically 
adjusted in order to maintain an approximately constant 
sensor signal /s- This makes our read-out signal virtu- 
ally unaffected by slow charge drifts and random charge 
rearrangements, which would otherwise prevent optimal 
bias. 

Fig. [Hb)(c) show schematic energy diagrams of elec- 
trons loading and unloading the quantum dot through 
the left barrier while the right barrier is raised high to 
completely cut off the channel. The electron tunnel rate 
is independently controlled by the left barrier gate volt- 
age Flb- In (b), both the high and low phase of the pulse 
are below the Fermi level Ey of the electron reservoir. 
Hence, there is no change in electron occupancy of the 
dot and the pulse lock-in detection signal /pulse is small, 
due only to capacitive coupling between the plunger gate 
and the sensor. Conversely, in (c), an electron can tunnel 
into the quantum dot from the reservoir during the high 
phase of the pulse (low in potential) and tunnel off when 
the pulse phase is low (high in potential). This change 



(a)- 



\ \ \ \ \ ^ ^ X \ v \ ^ 




< 20 



G3 








\ 




'^E3 1 






20 mV 



> 
E 



1.25 



»/p(V) 



1.28 1.25 







A/=2 \m' 





[arb. units) 

5 



FIG. 2: (Color online) (a) SET sensor lock-in current /s vs. 
Vp and Vrb, mapping out the charge stability plot of the 
quantum dot. Dynamic compensation is applied to the SET 
sensor, (b) Measured pulse lock-in signal /pulse from SET 
sensor with a symmetric square wave of peak-to-peak voltage 
V^puise = 2 mV applied to dot plunger P at /pulse =487 Hz. The 
first charge transition is not visible because the tunnel rate 
into/out of dot is too low with respect to /pulse- (c) Lock- 
in detection signal, with V^puise = 20 mV, at /puise=444 Hz 
for N — 2 ^ transition, extracted along the yellow dashed 
line in (d) at Vpuise = 20 mV. The red arrow indicates the 
orbital excited state, (d) Derivative of the detection signal 
with respect to Vp. Orbital ground and excited states are 
observed at the loading edge as white parallel lines. 



of electron occupancy in the dot induces a much larger 
lock- in detection signal /pulse- 



III. QD OCCUPANCY 

Fig. [2l^a) shows the measured QD charge stability 
map as a function of right barrier gate voltage Frb and 
plunger gate voltage Vp. As we reduce Vp^ the number 
of electrons in the quantum dot is reduced one by one 
until we observe no more charge transitions in the sta- 
bility map. The charging energy of the dot is around 
10 meV for high occupancy and increases significantly 
(up to 20 meV) in the few-electron regime. This indi- 
cates that the dot size is dramatically affected by the 
number of electrons at low N and strongly suggests that 
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FIG. 3: (Color online) (a) Lock-in signal from sensor S with 
V^puise=40 mV at /pulse =487 Hz applied to gate P, as a func- 
tion of Vp and Vlb- Here we apply a linear scaling of the 
lock-in output to directly indicate the excess electron occu- 
pancy AA^ on the dot. INSET: different region focusing on 
the first transition shown to enhance visibility of both excited 
state and unloading edge, (b) Orbital excited-state energy of 
the first 4 electrons at yLB=0.3 V.(c) Positions of AN = 0.5 
for ground states (GS) and first excited states (OS). These 
points are extracted from panel (a) where the tunnel rate 
from the reservoir to the dot is around 974 Hz. 



we have achieved the few-electron regime. In light of 
this, we believe that the last transition observed is likely 
to be ascribed to the last electron (A^ = 0). Simultane- 
ously, we measure the pulse lock-in detection signal /pulse 
when a symmetric square wave with Vpuise of 2 mV and 
/pulse =487 Hz is applied to gate P, and the resulting sta- 
bility map is plotted in Fig. [2fb). Identical charge tran- 
sitions are detected using /pulse , but with an improved 
signal-to-noise ratio compared with the sensor signal Is 
plotted in Fig. Efa). This allows us to detect a num- 
ber of extra transitions which are believed to originate 
from charge impurities nearby the sensor. Indeed, these 
features have different slopes with respect to the main 
ones and their signal strength is significantly reduced. 
This indicates that the originating charge displacement 
is elsewhere located in the substrate with respect to the 
gate-defined QD. 

Towards the one-electron limit, as the dot potential 
well becomes shallower and the tunnel barrier widens, the 
tunnel rate decreases. Here, the tunnel rate for the first 
electron falls significantly below 974 Hz (2x/puise Hz) 
and so the A^=0^1 transition is not visible in Fig. [2fb). 
Although electron tunneling is allowed energetically (see 
Fig.[2fa)), the high pulsing frequency used does not pro- 
vide an electron with sufficient time to tunnel onto/off 
the dot. In this case, the tunneling time is longer than 
1 ms. 

In order to observe the excited states of the dot, we 



increase Impulse to 20 mV as shown in Fig. [2]^c-d) for the 
= 2 ^ 3 transition. The pulse lock-in signal /pulse 
increases as soon as the three-electron ground state is 
pulsed below (green arrow in Fig. [2](c)). The sig- 
nal then decays slowly with increasing Vp indicating a 
decrease in tunnel rate as the ground state moves away 
from resonance with the Fermi level. When the three- 
electron excited state is pulsed below Ep the tunnel rate 
rises rapidly again, resulting in a second increase in the 
detection signal (red arrow in Fig. Efc)). When we in- 
crease Vp further, so A/" = 3 at both pulse levels, the 
signal falls low again since the tunneling events are ener- 
getically forbidden due to Coulomb blockade. Fig. [2fd) 
shows the derivative dIpxi\se/dVp of the pulse lock-in de- 
tection signal as a function of Vpuise and Vp. The left 
white line represents the ground-state loading edge. The 
second white line in parallel with the ground state is the 
excited state loading edge, while the unloading edge ap- 
pears as a black line. Here by loading (unloading) edge 
we refer to the potential configuration for which an elec- 
tron can start to occupy an energy level in the upper 
phase of the pulse (not) being allowed to empty it during 
the lower phase. 



IV. ORBITAL SPECTRUM 

We next study in detail the orbital excited states for 
the first few electrons by increasing Vpuise to 40 mV and 
plotting the charge stability diagram of the dot as a func- 
tion of Vi,B and Vp in Fig.[3ja). We identify these excited 
states as orbital states, as opposed to valley states, be- 
cause their energies are much larger than the valley split- 
ting in our device (~230 /ieV), as we shall show in the 
next section. In Fig. [31(a), we have converted the pulse 
lock -in signal /pulse 

to a corresponding average change of 
electron occupancy AN in the dot using a simple linear 
map. We identify AN = 1 as the level where /pulse satu- 
rates, which occurs when the tunnel rate significantly ex- 
ceeds the pulse frequency /pulse, as seen in Fig.[2](c) at Vp 
beyond the E3 riser. Note that AN = 2 where two tran- 
sitions overlap, as the pulsing level exceeds the charging 
energy allowing occupancy to change by two electrons. 
The ground states and the orbital excited states of the 
first four electron transitions are highlighted by green 
and red lines, respectively. We also observe a decrease 
in tunnel rate as we decrease Flb (thus increasing the 
barrier height), resulting in AA^ ^ as the tunnel rate 
falls below /pulse- In Fig- ^^)^ we plot the first orbital 
excited state energy £^orb with respect to the electron 
number A^, extracted from Fig. [31(a). Note that the en- 
ergy conversion factor, a ~ = 0.3 eV/V, is derived from 
considerations relevant to valley splitting measurements, 
as discussed next. As Vp is reduced, the dot becomes 
smaller leading to an increase in the orbital level spacing 
from 1 meV for A/" = 4, to 8 meV for A/" = 1. 

In Fig.[3](c), we plot the positions in gate voltage space 
(^lb; ) corresponding to AA^ = 0.5 for the ground and 
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orbital excited states = 1 to 5. These are the points 
where the tunnel rates are comparable to the pulse fre- 
quency /pulse- As the dot occupancy increases, a greater 
barrier height (or lower Vi,b) is required to maintain a 
constant tunnel rate. We note that both the ground 
states (GS) and the orbital excited states (OS) follow 
regular trend lines, except for the N = 5 ground state, 
which appears more like an orbital excited state. This is 
due to the ground orbital levels being full with electrons 
and starting to fill the next orbital. 



V. SPIN- VALLEY INTERACTION 

We now study the spin filling of the first four electrons 
into the valley states of the quantum dot by examining 
the Zeeman shifts of the ground and excited states in a 
magnetic field parallel to the oxide interface in the range 
-6 T < 5 < 6 T, as shown in Fig. IH^a). Here we set 
Vpuise = 20 mV and plot the derivative of the pulse lock-in 
current, dlp^ise/dVp. The data reveals a number of levels, 
including spin ground and excited states, for what appear 
to be two distinct valley levels, as we explain below. 

We first consider the unloading edge of the charge tran- 
sitions, which appear as dark lines in Fig. IH^a). For the 
first transition {N = 0^1), the unloading edge moves 
towards less positive Vp with increasing magnetic field 
indicating that the first electron is spin-down ||). 
Assuming that the g-factor = 2 for electrons in silicon 
and using the Zeeman energy equation ^g/j^sB^ where 
/iB is the Bohr magneton, we fit the energy change of 
the unloading edge. This results in a plunger gate volt- 
age to energy conversion factor, a ~ 0.3 eV/V. This is 
consistent with values for a similar quantum dot reported 
in Ref.^^. Note that the evaluation of the a-factor from 
the Zeeman splitting is a robust approach in Si-based 
systems. Because of the weak spin-orbit coupling^^, only 
very small deviations (less than 1%) of the g-factor from 
the free-electron value are observed in Si, as opposed to 
the case of III-V semiconductors where much larger de- 
viations are found. Next, we analyze the unloading edge 
of the second transition (TV = 1^2). The energy level 
first increases with magnetic field for \B\ < 2 T, indi- 
cating that a spin-up electron was unloaded, and then 
decreases with \B\ above ~ 2 T. This indicates that the 
N = 2 ground state changes from a singlet -i=(|H) — 

at low 5, to a triplet |H) for |5| > 2 T. A similar kink at 
\B\ ^ 2 T (marked by a yellow arrow) is observed for the 
third transition (A^ = 2 ^ 3), although here the unload- 
ing edge moves downwards at low \B\ and then upwards 
at high \B\. Finally, the unloading edge of the fourth 
transition {N = 3^4) rises linearly with consistent 
with a spin-up electron unloading for all magnetic fields 
up to 6 T. 

In order to explain these data, we use a model based 
on two non-degenerate valley states (Vi, V2) whose spin- 
states depend on B field according to the Zeeman en- 
ergy equation (Fig. IHb)). This model relies on the fact 
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FIG. 4: (Color online) (a) Grey-scale plot of the QD differen- 
tial occupancy (A'^N) in magnetic field for the first 4 electron 
transitions at Vpuise =20 mV and /puise=487 Hz. Orbital ex- 
cited states, Zeeman splittings and valley-orbit splitting are 
observed. Vp represents arbitrarily shifted Vp. (b) Schematic 
of the evolution of two non-degenerate valley eigenstates for 
positive B-field. Vi (Vi') and V2 (V2') represent loading 
(unloading) edge states. Each valley state Zeeman splits into 
two levels at finite B-fields. The blue dots represent already- 
occupied electron states. 



that bulk silicon has six degenerate valleys, but when 
confined to two dimensions, such as a Si/Si02 interface, 
these separate into four valleys with high effective mass 
and two lower-energy F- valleys^. The F- valleys are also 
non-degenerate with a splitting AEy = Ey2 — E-yi that 
is typically well below 1 meV. In previous ground-state 
magneto-spectroscopy studies on a similar dot we found 
AE^J - 0.1 me via, which we note is much smaller than 
the first orbital excited state energies (^orb = 1—8 meV) 
as long as dot occupancy is kept low {N < 4). From 
Fig.m^b), we see that the first electron fills the lower val- 
ley Vi as spin-down for all B. When the second electron 
loads (A^ = 1 ^ 2 transition) it will fill Vi as spin-up at 
low values of forming a spin singlet --^(|H) — |ti)), 
however at large when the Zeeman energy exceeds 
the valley splitting A£^v, it will preferentially load into 
the upper valley V2 as spin-down, to form a triplet state 
1^^). Moving to the third electron, we see that it will fill 
as spin-down in valley V2 for low and then as spin- 
up in valley Vi once the Zeeman energy exceeds AEy. 
The fourth electron always fills as spin- up in valley V2, 
unless there is a nearby orbital state (within the Zeeman 
energy), which is likely to be the case for higher dot oc- 
cupancy than currently reported {N > 4). In Fig. HJ^a) 
we can clearly see the orbital excited states (white lines) 
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between the loading and unloading edges. They predom- 
inantly move downwards with the Zeeman energy, imply- 
ing that they load spin-down electrons. 

We find that the level crossings in Fig. |4] (marked by 
yellow arrows) occur at \B\ ~ 2 T for N = 1..3, indicat- 
ing that the valley splitting AEy = gfisB ~ 230 fieV 
remains approximately constant for small electron occu- 
pancy. This value is roughly double the valley splitting 
we observed in a previous device^^, most likely due to 
the larger electric fields employed in the present study 
The observation of a valley splitting that is independent 
of electron number confirms the prediction^ ^ that the val- 
ley exchange Coulomb integral is negligible. 

The ground-state loading edges (lower white fines in 
Fig. m^a)) for transitions 2 to 4 show the same trends in 
magnetic field as the unloading edges. However, for the 
first transition (TV = ^ 1) the ground state (marked by 
a green dashed line) is not visible. We know this state is 
present because the distance between the loading and un- 
loading edge must equal Vpuise = 20 mV. We believe this 
is due to a very low (loading) tunnel rate, which we pre- 
viously observed for low electron occupancy in Fig. [3l^a). 
Note that the observation of the spin excited state is com- 
patible with the pulse frequency in use which is chosen to 
be higher than the relaxation rate. Indeed, the relaxation 
rate in the window of magnetic field applied is expected 
to be in the range 0.1-10 Hz^^^^, whereas we use pulse 
sequences of few hundreds of Hz. This accounts for the 
kink at 2 T formed by the lower valley spin up state and 
the upper valley spin down state. Finally, we consider 
the = 5 ground state, which has a tunnel rate that 
follows the trend for the orbital excited states for N = 1 
to 4, as shown in Fig. [3jc). We can now understand this 



in terms of a shell structure, where the first four elec- 
trons fill the spin and valley states of the lowest orbital, 
and the fifth electron occupies the next available orbital 
level. This is consistent with a shell structure of = 4 
observed in other silicon quantum dot^ii^. 



VI. CONCLUSIONS 

In conclusion, we have presented excited state spec- 
troscopy of a nearly-closed silicon quantum dot using 
charge sensing and a pulsed-gating technique, thus en- 
abling clear identification of the spin, valley and orbital 
states for the first four electrons. As the occupancy in- 
creased from = 1 to 4 electrons, we found that the 
valley splitting for the lowest orbital level remained ap- 
proximately constant at 230 /ieV, while the next orbital 
level energy decreased from 8 meV to 1 meV. Given the 
increasing interest in quantum information processing us- 
ing spin and valley states in silicon quantum dots^^, it is 
important that the multi-valley level structure of these 
systems is well characterised experimentally. 
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